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ABSTRACT 

Context. Direct observations of accretion disks around high-mass young stellar objects would help to discriminate between different 
models of formation of massive stars. However, given the complexity of massive star forming regions, such studies are still limited in 
number. Additionally, there is still no general consensus on the molecular tracers to be used for such investigations. 
Aims. Because of its close distance and high luminosity, IRAS 17233-3606 is a potential good laboratory to search for traces of 
rotation in the inner gas around the protostar(s). Therefore, we selected the source for a detailed analysis of its molecular emission at 
230 GHz with the SMA. 

Methods. We systematically investigated the velocity fields of transitions in the SMA spectra which are not affected by overlap with 
other transitions, and searched for coherent velocity gradients to compare them to the distribution of outflows in the region. Beside 
CO emission we also used high-angular H2 images to trace the outflow motions driven by the IRAS 17233-3606 cluster. 
Results. We find linear velocity gradients in many transitions of the same molecular species and in several molecules. We report the 
first detection of HNCO in molecular outflows from massive YSOs. We discuss the CH 3 CN velocity gradient taking into account 
various scenarios: rotation, presence of multiple unresolved sources with different velocities, and outflow(s). Although other inter- 
pretations cannot be ruled out, we propose that the CH 3 CN emission might be affected by the outflows of the region. Higher angular 
observations are needed to discriminate between the different scenarios. 

Conclusions. The present observations, with the possible association of CH 3 CN with outflows in a few thousands AU around the 
YSOs' cluster, (i) question the choice of the tracer to probe rotating structures, and (ii) show the importance of the use of Ho images 
for detailed studies of kinematics. 
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QO ' "I ■ Introduction Over the last years, indirect evidence of the existence of ac- 
f^) , cretion disks have been provided by observations of phenom- 
■ ' 1 ena such as collimated outflows and jets with properties sim- 
g ■ Our understanding of the formation of massive stars has largely ilar t0 those originating from low-mass YSOs, thus implying 
^ ; improved in recent years. From a theoretical point of view, it has a similar f orma tion mechanism for low- and high-mass stars. 
^ . been demonstrated through 2-D and 3-D radiation hydrodv- However, direct observations of accretion disks in high-mass 
I namic simulations (|Krumholz et al| |2009t |Kuiperet al| |2010|), sta r-forming regions are still limited to few cases, mostly around 
m > , that radiation pressure is not a barrier to form stars with masses YSOs with masse s < 20 M Q and luminosities < 10 4 L (e.g., 
km ■ >8 M Q and spectral type O and B and that they form through iKraus et allfloToh. The situation is less clear for m ore luminous 
^ ; disk accretion. In fact, the two main theoretical scenarios pro- object s (ICesaroni et all 120071) . although recently IZaoata et all 
03 ■ posed to explain the formation of massive star s are accretion- (pJrM |201 Obi) suggested the presence of a Keplerian infalling 
based mech anisms: (1) the core accretion model (|McKee & Tatj ring around a central ob j ect of at least 60 Mq . This casts some 
|2002l [2003, where a massive star forms from a massive core dou bts on the existence of disks around O-type stars as no con- 
fragmented from the natal molecular cloud and (2) the com- clusive evidence for coherent, stable accretion disks around them 
petitive accretion model dBonnell et al| |2007|), where a molec- as been preS ented so far. The clustered mode of high-mass star 
ular cloud initially fragments into low-mass cores, which form formation and the typical large distances of high-mass star form- 
stars that compete to accrete mass from a common gas reservoir. ing reg ions challenge direct observation of accretion disks. In 
Both models predict the existence of protostellar accretion disks addition, it is still unclear which molecular lines are best suited 
around massive young stellar objects (YSOs), and therefore the for suc h studies. They should have low optical depths, and trace 
presence of molecular outflows, although the properties of both nly the hot gas around the central protostar, and not its enve- 
disks and outflows could be different. In fact, the competitive ac- lope . In the past, molecules such as NH 3 , CH3CN, HCOOCH3, 
cretion model suggests that massive stars should always form in and C 34 S have been p r oposed for this purpose (i.e IZhang et all 
densely clustered environments, and so the circumstellar disks [T^. ICesaroni et all fl997[ fl999l iBeuther et all l2004 l2009i. 
could be severely perturbed by interaction with stellar compan- 
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However, there is still no general cons ensus on the molecular 
tracer to use (i.e., iBeuther et al.l 1200 5b). 

IRAS 17233-3606 (hereafter IRAS 17233) is one of the 
best lab oratories to s t udy th e problem of massive star for- 
mation. Faun dez et al.l d2004l) reported a luminosity of 1.7 x 
10 4 L Q obtained by integrating the spectral energy distri- 
bution and assuming d = 0.8 kpc. Distances between 
700 pc and 2.2 kpc are reported in the literature (i.e., 
iMiettinen etaE 120061: iForster & Caswelli 119891) . However, the 
latter is based on maser emission lines . Using the rotation 
curve determined by Brand & Blitz ( 1993) and velocities based 
on thermal line s (Bronfman et al., 1996; M acLeod et al l fl998t 
Miett inen et all 120061: iLeurini et all 1201 ll) . the kinematical dis- 
tance of IRAS 17233 is constrained to < 1 kpc. The far distance 
(~ 15 kpc) is unlikely because of the source extremely bright 
continuum and line emission at practically all observed wave- 
lengths, which would indicate exceedingly high luminosities if 
the source were at the far distance. In the following, we assume 
a distance of 1 kpc for IRAS 17233. 

Thanks to co mbined high angular resolution observations of 
different tracers, ILeurini et al.l (120091) found evidence for mul- 
tiple outflows driven by different hyper-compact H n (HCH n ) 
region s resolved at milliarcsec-scale with VLA by Zapat a et al.l 

(2008) (see Fig. [TJ. Their parameters are typical of massive 
YSOs (early B-type ZAMS stars) and in agree ment with the 
measured bolometric luminosity of the source. IBeuther et al.l 

(2009) found evidence for a large-scale rotating structure de- 
tected in high excitation lines of NH3 along the ESE-WNW di- 
rection. On th e other hand, t he hi gh spatial-resolution maps in 
CO and H 2 by ILeurini et' al. (2009) indicated a compact bipolar 
outflow located along the NW-SE, thus questioning the interpre- 
tation of the NH3 gradient as a rotating structure. 

We present here observations of the massive star forming re- 
gion IRAS 17233 in several molecular tracers around 230 GHz 
taken with the Submillimeter Array. The idea of this study is 
to analyse the spatial distribution and velocity field of different 
lines, compare them to the distribution of outflows in the region, 
and verify whether we find evidence for rotating structures per- 
pendicular to the outflows or not. Given the large bandwidth of 
SMA spectra (~ 4 GHz) and the richness of the molecular spec- 
trum of the source, such observations allow us to study the dis- 
tribution of several lines, and not only of the "best" disk tracer 
candidates. 

2. Molecular outflows from the YSOs in 
IRAS 17233-3606 

In order to investigate the presence of possible accretion disks 
around the protostars in the inner region of IRAS 17233, a de- 
ta iled knowledge of th e molecular outflows in the area is needed. 
In lLeurini et al.l d2009l) we analysed the emission of CO isotopo- 
logues observed simultaneously to the data presented in this pa- 
per (see Sect. [3j> and compared their distribution to high-spatial 
resolution H2 observations. The main results of that study can be 
summarised in the following points (see Fig. 1): 

1 . The data reveal a clumpy extended structure with well sepa- 
rated blue- and red-shifted emission, and an overall structure 
roughly aligned along the N-S direction. The outflow is as- 
sociated with extremely high velocities, up to 200 and 

+ 120 km s _1 with respect to the ambient LSR velocity. 

2. Using high angular resolution maps of H2, cm-continuum, 
and maser emission, multiple outflows can be distinguished 
in the region. One (called here OF1) is compact (~ 5"-10") 



Table 1. Observational parameters. 
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and located along the NW-SE direction, as detected in CO 
and SO: it shows extremely high-velocity red-shifted emis- 
sion and it is associated with a jet traced by H2 emission as 
well as with red and blue H2O maser spots. The extended 
emission along the N-S direction is probably due to multi- 
ple outflows. The picture derived by connecting the position 
of the CO high velocity clumps suggests two possible main 
axes, here labelled as OF2 and OF3 and drawn in Fig. [T|by 
dashed lines. Interestingly, OF2 is associated with a coun- 
terpart on smaller scales (~ 5", corresponding to ~5000 AU 
at D — 1 kpc) traced by red- and blue-shifted OH masers 
(Fig. 1 1(b)) . In addition, another possible outflow direction is 
indicated in Fig.Q]by the dotted line. 
3. The powering sources of the flows remain unidentified, 
as they are centred on a region where 4 HCH n regions 
(VLA 2a, b, c and d, see Fig. |l(b)| and |2]i have been observed 
on a sub-arcsec scale. 

3. Observations 

We observed IRAS 17233 with the SMA interferometer on April 
10, 2007 in the compact configuration with seven antennas. The 
receivers operated in a double-sideband mode with an IF band of 
4-6 GHz so that the upper and lower sideband (USB and LSB, 
respectively) were separated by 10 GHz. The central frequencies 
of the upper and lower sideband were 220.4 and 230.4 GHz, 
respectively. 

The observations presented in this stud y were performed si - 
multaneously to those of CO discussed by ILeurini et al.l J2009). 
For details on the weather conditions, calibration, and reduction 
of the dataset we refer to their Sect. 2. We checked the absolute 
positional uncertainty on the quasar 3C273, our phase calibrator, 
and found it to be < 078 for the right ascension and < 0V5 for 
the declination. 

The resulting synthesised beams, the velocity resolution used 
in the analysis and the typical r.m.s of the maps are given in 
Table Q] The centre of the maps is at a 2 ooo=17' ! 26 ffl 41' ! .757, 
<520oo=-36°09'0'.'50. The conversion factor from flux density 
to brightness temperature in the synthesised beam is ~2.25 
K (Jy/beairfT 1 for the lower side band, and ~2.53 K (Jy/bearrfT 1 
for the upper side band data. 

4. Observational results 

4.1. Continuum emission 

The continuum emission of the region at 1.3 mm is shown in 
Fig. [2] This was derived using only the upper side band data, 
which have a lower number of lines than the lower side band 
data and therefore line subtraction is less problematic. The peak 
of the emission is found close to the cluster of four co mpact 
sources detected at cm-wavelengths by Zap ata et al.l (l2008). This 
also corresponds to the so-called "maser zone", where strong 
H2O, CH3OH, and OH masers a re det ected (e.g. JCaswell et all 
1980; iFixetall 119821: iMentenl 1199 ll see Fig. [T(b)l ). Weak (~ 
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Fig. 1. Left panel: Integrated emission of the blue- and red-shifted wings in the CO(2-l) line overlaid on the H2 emission at 2. 12/vm 
(grey scale). The solid contours show the EHV blue- (v=[-200,-130] km s -1 ) and red-shifted emission (v=[90,120] km s -1 ); the 
dashed contours mark the HV blue- (v=[-130,-25] km s -1 ) and red-shifted emission (v=[16,50] km s -1 ). The dashed lines and the 
solid arrow outline the possible molecular outflows identified in CO (OF1, OF2, OF3, see discussion in Sect. 0); the dotted line 
marks the alternative direction of the OF2 outflow. The central box outlines the region mapped in the right panel. The squares 
mark the positions of the red- and blue-shifted CO(2-l) peaks. Right panel: Continuum emission at 1.3 mm from the maser zone in 
IRAS 17233 (solid black cont ours from 50% of peak emission in step of 1 5%, see S ect.|4JJi. Overlaid on the map are the positions of 
the OH (blue and re d squares. iFish et al.ll2005l) . CH3OH (black triangles. IWalsh et"aQ [l998) and water (blue and red open triangles, 
IZapata etal.l l2008) masers. The gree n contour levels show the continuum emission at 1.3 cm (from 1.4 mJy beam ', 3cr, in step of 



1.4 mJy beam Zapata et~aQ 120081) . while the grey scale map is the H2 emission. The solid blue and red contours represent the 



SO(56 - 45) blue- and red-shifted emission (blue: from 80 to 104 Jy km s -1 beam -1 in step of 10 Jy km s -1 beam -1 ; red: 10 and 
12 Jy kms -1 beam -1 ). The dashed red contour shows the EHV red-shifted CO(2 - 1) emission. OF1, OF2 and OF3 are marked as 
in left panel. 



3cr) emission is also detected from the H 11 region on the east of 
the maser zo ne, already mapped at cm-wavelengths by several 
authors (e.g jHughes & MacLeod! 1 1 9931: IWalsh etal.Ul998l) . 

We performed a 2-D Gaussian fit of the 1.3 mm con- 
tinuum emission, and determined its peak to be at a^ooo = 
\l h 26 m A2A55\ 6 nm = -36°09'18.047" with positional errors 
of 0708. The peak flux (F v ) is ~ 2.09 + 0.07 Jy beam -1 , the 
integrated flux (S v ) 7.20 Jy. However, given the richness of the 
molecular spectrum of the region at this frequency, the contin- 
uum emission may have been overestimated because of line con- 
tamination. The inferred source size is 67 1 + 072 x 57 1 + 072 
with a P.A. of 64° .2 + 7°.5, corresponding to a deconvolved size 
of 573 x 277, and a P.A. of -77°.2 (see TableEJ. 

The peak of the mm continuum emission do es not corre- 
spond to any of the compact sources detected by Zar>at a~"et al.l 
(2008) at 1.3 cm, but it is shifted to the west. VLA 2a, b, c and d 
are all within ~ 2" from the peak of the mm continuum emission. 
However, the displacement between the mm continuum peak and 



the strongest cm source VLA 2a is only ~ (076, 070), while the 
offset from the CH3OH mas er spots (whose u ncertainty in the 
absolute position is ~ 1". IWalsh et"aU 1 1998b is ~ (174,073). 
Given the uncertainty in the absolute position of our SMA data 
(^ 078), this may not be significative. 

We can calculate beam-average d gas column density and 
mass using the equatio ns given by Beutheretal. (2002) cor- 
rected in the erratum by [Beuther et al.l (12005 al) . We assume the 
authors' default values for the grain size (0.1 pm), grain mass 
density (3 g cm -3 ), gas-to-dust ratio (100) and a typical dust 
temperatures of 200 K for hot cores. For the grain emissivity, 
we used Q v = 7.5 x 10~ 4 (125pm/A'f,/3 = 1.6, as given by 
Beut her et al.l (12005 ah . We derived a beam-averaged H2 column 
density of l.OxlO 24 cm -2 , and a mass of 12 M Q . However, the 
mass and column density derived in this fashion are affected by 
the uncertainties in the peak and integrated flux determinations 
due to line contamination. Based on unpublished SMA observa- 
tions at 217 GHz in the same configuration of the array as for 
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Fig. 2. Continuum emission at 1.3 mm (black contours) and 
1.3 cm (red contours) of IRAS 17233. Contours are from 
2 mJy beam -1 (5cr) in step of 2 mJy beam -1 (1.3 cm emis- 
sion), and from 0.12 Jy beam -1 (3<x) in steps of 0.24 Jy beam" 1 
(1.3 m m emission). The fo ur HCHn VLA 2a, b, c and d identi- 
fied bv lZapataetatl d2008l) are also labelled. 



the current data (Leurini et al. in prep.), we derived a second in- 
dependent measurement of S v and F v at a wavelength close to 
that presented here. The peak flux at 217 GHz is estimated to 
be 2.10 + 0.04 Jy beam -1 (on a beam of 3'.'31 x 2'.'54), while 
the integrated intensity of the source is 6.69 Jy. These values are 
very close to those derived at 230 GHz, with an uncertainty of 
less than 10% in the peak flux (but on a different beam) and on 
the integrated intensity. Therefore, the uncertainties in the esti- 
mates of the mass and of the column density are dominated by 
the uncertainty in the dust temperature. 

4.2. Line emission 

The molecular spectrum associated with the peak of the 1.3 mm 
continuum emission is very rich in lines, and typical of hot cores 
in massive YSOs. The LSB spectrum of the source (Fig. [3) is 
dominated by the CH3CN (12^ - 11k) band: transitions up to 
K — 9 are detected, together with lines from vibrationally ex- 
cited levels (up to K — 7), and from the CH^ 3 CN isotopologue 
(up to K — 6). This points to high excitation conditions since 
the lower level energies of these lines range up to 840 K. The 
detection frequency of lines in the USB spectrum is less than 
for the LSB band, and the band is dominated by the CO (2-1) 
transition. 

Given the complexity of the region and the richness of the 
spectrum, line identification in molecular spectra of hot cores is 
challenging. In the analysis of the spectrum of IRAS 17233, we 
made use of the XCLAS S prograrrQ. The me thod used in the 
program is described by IComito et al.l (120051) (and references 
therein), and it is based on a simultaneous fit of all lines in a 
spectrum under the LTE approximation. The model takes into 
account the beam dilution, the line opacity, and the line blend- 
ing, and produces a synthetic spectrum of the source for given 
physical parameters (7\j n , N spec i es , source size). The molecular 
spectroscopic parameters are taken from the CDMS catalogue 
(Mtil ler et aUl200ill2005l) and from the molecular spectroscopic 
datab ase of the Jet Propulsion Laboratory (JPL, see lPickett et al.L 
Il998h . 



1 http://www.astro.uni-koeln.de/projects/schilke/XCLASS 



Although this approach delivers physical quantities such us 
rotational temperature, column density and source size for each 
of the molecular species included in the model, we believe that 
the complexity of the region is too high to infer meaningful re- 
sults with such a simple approach. The analysis of the hot core 
chemistry of the source is beyond the scope of this paper and will 
be treated in a dedicated study. Nevertheless, a LTE synthetic 
spectrum, obtained with plausible input values for the physics of 
the gas of the inner region around the protostars, is fundamental 
for a correct identification of spectral features, since it predicts 
intensities for any transition of a given molecular species at all 
known frequencies. Based on the line profiles, and on the distri- 
bution of several transitions which show extended emission (see 
Sect. 15. II ), we assumed two components inside the beam, the first 
with a compact (< 2") source with a temperature of ~200 K and 
a line width of 4 km s -1 , and a second more extended (~ 10") 
source with lower temperature (50-70 K) and broader line width 
(~ 10 km s -1 ). Table|2]lists all transitions identified in the LSB 
and USB spectra of IRAS 17233 towards the peak of the mm 
continuum emission, which are analysed in this paper, and their 
line parameters obtained by using Gaussian profiles. 

5. Data analysis 

The frequency setup of our observations covers the CH3CN 
(12k - Hk) band, used in the past to trace disks and/or toroids 
around massive YSOs, as well as species which have been de- 
tected in molecular outflows from low- and high-mass protostars 
(i.e., H 13 CO, CH3OH, OCS). In addition, several high energy 
transitions are detected, which are potentially good candidates 
to search for tracers of rotation. 

To study the velocity field traced by different molecular fea- 
tures, we derived maps of the zero-th and first moments of dif- 
ferent transitions. We limited our analysis to lines which are not 
affected by overlap with other transitions. This was verified with 
the help of our LTE model (see § 14.21) . However, since several 
spectral features remain unidentified, it is possible that lines used 
to derive the distribution of the zero-th and first moments are af- 
fected by emission from species not included in our model. 

Table [2] lists all transitions identified in the LSB and USB 
spectra of IRAS 17233 towards the peak of the mm continuum 
emission and used to derive the zero-th and first moment maps. 
Part of the LSB spectrum towards the peak of the continuum 
is shown in Fig. [3] as example of the complexity of the source. 
Figures l4lfT0l and [T3lfT4l show the derived distributions of the 
zero-th order and first moments computed over a region where 
the signal-to-noise ratio in the line data cubes is larger than 5cr 
(corresponding to 0.5 Jy beam -1 ). 

5.1. Spatial distribution 

As a first step to investigate the association of given transitions 
with the molecular outflows in the region, with potential accre- 
tion disks or with the envelope around protostars, we derived the 
integrated intensity of different lines to study their spatial distri- 
bution. In Figs.|U|5]|6] and[7]we present the integrated intensity 
distribution of selected transitions (Table All the maps are 
elongated along the PA. of the beam of the observations. For 
all maps, we performed a 2-D Gaussian fit to infer the FWHM 
size and peak position of the integrated intensity distributions. 
The results are listed in Table [3] All emission maps are resolved 
with our current resolution, except those of the HNCO (IO3-93), 
HC3N (24_ 2 - 23 2 ) v 7 = 2, CH 13 CN (12 2 - 11 2 ), and CH3CN 
(12 2 — 11 2 ) Vs = 1 lines. The fit suggests that the peak of the 
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Fig. 3. Lower side band continuum-subtracted representative spectra of IRAS 17233 at the peak of the 1 .3 mm continuum emission. 
All transitions used to derive the zero-th a nd first moment ma ps of Figs. l4lfT4l are labelled. For reference, we also mark the C ls O 
(2 - 1) and SO(5& - 5s) lines (discussed in iLeurini et all 120091) . and the ladder of the CH3CN (labelled with the corresponding K 
quantic number) and CHi 3 CN (12-11) lines (labelled with the corresponding K quantic number followed by 13 C). 



molecular emission does not coincide with the 1 .3 mm contin- 
uum maximum, but it is located slightly to its NE, with a mean 
offset of (9'.'0, -1773) from the centre of the maps, (fX'6, (X'3) 
from the peak of the mm continuum emission. This could indi- 
cate that the hot core i s associated with th e compact cm source 
VLA-2a identified by IZapata et all (120081) (also shifted to the 
east of the mm continuum), whereas the bulk of continuum is 
tracing colder material. Opacity effects could play a role in this 
displacements. Actually, the CH3CN spectral pattern shows sim- 



ilar brightness temperatures for lines with K < 5 (Fig. [3}, sug- 
gesting non negligible opacities. However, the distribution of 
less abundant isotopologues (see 13 CS) or vibrationally excited 
transitions (HC3N V7 = 2 lines) seem to peak closer to the con- 
tinuum peak, but still shifted to the east, near VLA-2a. 

Methanol shows extended emission along the axis of the 
OF2-OF3 outflows, and an elongated structure towards NW, i.e. 
in the direction of the OF1 red-shifted outflow detected in H2, 
extremely high velocity CO and SO (see Fig. |1 (a)) > and indicated 
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Table 2. Spectral features detected in the IRAS 17233 spectrum and analysed in the present paper 
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CH 3 CN(12 2 - 11 2 ) v 8 = 1 


221367.67 


596 


-3.63 ± 


0.37 


8.14 ±0.97 


178.5 ± 17.3 


CH3OH-E (8_i-7 ) 


229758.76 


89 


-3.22 ± 


0.02 


7.94 ± 0.04 


484.0 ±2.1 


13 CS(19- 18) 


230317.53 


99 


-2.50 ± 


0.04 


6.36 ±0.11 


118.9 ± 1.7 


OCS (19-18) 


231060.99 


100 


-3.04 ± 


0.02 


10.51 ±0.04 


574.2 ± 1.7 


CH 3 OH-A (10 2 - 9 3 ) 


231281.10 


165 


-3.72 ± 0.03 


7.31 ±0.08 


165.6 ± 1.3 



Table 3. Result of the 2-D Gaussian fit of the integrated intensity distributions. 



transition 


Aa a 


A<5 a 


FWHM observed size b 


FWHM deconvolved size 


cont. 


8'.'455 


-177 551 


67 1 x 5 


'1 


573 x 277 


C 2 H 5 CN(24 22 , - 23 221 ) 


8'.' 930 


-177 5 38 


579 x 3 


'2 


276 x 272 


HNCO(10 3 - 9 3 ) 


8'.' 987 


-177140 


570 x 1 


'9 




HC 3 N (24 - 23 ) v 7 = 2 


8'.' 941 


-177421 


575 x 2 


'3 


175 x 074 


HC 3 N (24_ 2 - 23 2 ) v 7 = 2 


8'.' 832 


-177475 


574 x 2 


'2 




HNCO(10o-9 ) 


8'.' 830 


-177 399 


579 x 2 


'9 


275 x 271 


CH 3 OH-£ (8 -7_i) 


97 091 


-177108 


579 x 3 


'9 


376 x 272 


CH 3 CN(12 8 -11 8 ) 


8'.' 885 


-177245 


575 x 2 


'6 


270 x a/4 


CH 3 CN(12 7 - 11 7 ) 


97010 


-177 305 


575 x 2 


'9 


272 x 170 


CH> 3 CN(12 2 - 11 2 ) 


8'.' 900 


-177400 


574 x 3 


'0 




C 2 H 5 CN (25 2 , 24 - 24 2 , 23 ) 


97076 


-177459 


578 x 3 


'1 


276 x 270 


CH 3 CN(12 4 - 11 4 ) 


97 113 


-177 3 86 


579 x 3 


'8 


373 x 272 


CH 3 CN(12, - ll 2 )v 8 = 1 


87975 


-177 375 


573 x 2 


'4 




CH 3 OH-£ (8-1 - 7 ) 


9! 247 


-167994 


578 x4 


'7 


474 x 279 


13 CS(19- 18) 


87850 


-177 559 


572 x 3 


'6 


371 x 177 


OCS (19- 18) 


97216 


-177424 


577 x 4 


'4 


470 x 279 


CH 3 OH-A(10 2 -9 3 ) 


97173 


-177 082 


573 x 3 


'8 


373 x 270 



a offset from the centre of the maps at a 200 o=17''26'"41 s .757, <5 20 oo=-36 09'0750 (see Sect. 3). 
b typical errors on the fit are of the order of 072. 



by an arrow in Figs.|4]|7] Interestingly, the occurrence of CH3OH 
emission along the OF1 axis is more evident in the two lower ex- 
citation transitions (reported in Table [2). As shown by the spec- 
tra in Fig. [3] this is not due to low dynamical range, since all 
methanol lines used in this analysis have a signal-to-noise ratio 
larger than 40 at the peak of the continuum emission. 

Also the OCS emission, and to some extent 13 CS emission, 
is elongated in the direction the OF2-OF3 outflows, and shows 
a wisp of emission in the direction of OF1 (Fig. [5}. 

Figure [6] shows the zero-th moment maps derived from dif- 
ferent lines of CH3CN. The spatial distribution of the 12 4 - 1 14 
line, the lowest CH3CN excitation line used in our analysis, ex- 
tends in the direction of the OF1 jet and along the axis of the 
OF2-OF3 outflows, while the other transitions are more com- 
pact. Also in this case we can exclude a bias due to different 
signal-to-noise ratios for the different CH3CN lines. 

The C2H5CN and HNCO (10 - 9 ) line maps are also elon- 
gated along the direction of the OF2-OF3 axes, but their emis- 
sion is more compact than for the other species discussed above. 



In conclusion, the spatial distribution of CH3OH, OCS, and 
CH3CN emission could be affected by the presence of the out- 
flows. In particular, for CH3OH and CH3CN we find, besides 
an elongation along the OF2-OF3 axes, also a good correlation 
with the direction of the OF1 outflow. Emission at the position 
Rl, the peak of the CO (2 - 1) red-shifted emission at extremely 
high velocity in OF1, is detected from the CH3OH (8_i -7q) and 
(8 -7i) lines and the CH3CN (12 4 - 1 1 4 ) transition. In Sect. EH 
we will analyse the first moment maps to investigate the kine- 
matics of the gas and better constrain the association of molec- 
ular emission with the molecular outflows from the IRAS 17233 
cluster. 



5.2. Velocity distribution 

For transitions analysed in Sect. 15. II with emission larger than 3 
beam widths along the minor axis of the beam, we also computed 
the first moment distribution. This excludes the HC3N transi- 
tions from the analysis. 
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The detected velocity shifts are of the order of Av ~ 
6 km s , i.e. values larger than velocity resolution of the data 
cubes (0.5 km s -1 ). All CH3CN maps are blue-s hifted compared 
to the systemic velocity (V sys ) of -3.4 km s -1 (Bronfman et all 
[l996h . However, there is a large uncertainty in the determination 
of V S ys (between -3.7 and -2.5 km s _1 from single dish N2H + 
and C n O observations. iLeurini et al.lfeoi ll) . probably due to the 
outflow multiplicity in the region. 

As suggested by the analysis of the integrated intensity dis- 
tributions, CH3OH clearly presents a velocity gradient along 
the direction of the OF2-OF3 outflows, with blue-shifted emis- 
sion to the north-east and red-shifted emission to the south- 
west. Additionally, the CH 3 OH (8_i - 7 ) and (8 - 7\)-E lines 
show red-shifted emission associated with the red lobe of the 
OF1 jet (see Fig. |l(b)) , A wisp of red-shifted emission towards 
OF1 is detected also in the CH3OH (10_ 2 - 9_ 3 ) line although 
marginally. 

For simplicity, we present the first moment distribution of the 
CH3CN (12 4 - 11 4 ) line in Fig. while the CH 3 CN (12 7 - 11 7 ) 
transition is presented in Fig. [10] We see the same behaviour 
seen in the CH3OH (10_2 -9_3) line in the (H4- 1 14) transition, 
the lowest energy analysed line of methyl cyanide: a linear ve- 
locity gradient along the E-W axis plus a red-shifted wisp of gas 
extending to the north-west towards the red-shifted lobe of the 
OF1 outflow. The CH3CN (12j - II7) line has a more compact 
emission, and shows a linear velocity gradient along the E-W 
direction, almost perpendicular to the OF2-OF3 axis. Although 
not shown here, the velocity field of the CH3CN (12g - lis) line 
agrees with that of the CH3CN (I27 - 1 17) transition. 

In order to further investigate the association of the CH3CN 
(I24 - II4) line with the OF1 outflow, in Fig.|9]we compare its 
first moment distribution to the H2 distribution and to other trac- 
ers of kinematics in the region. The SO blue- and red-shifted 
emission and the blue- and red-shifted H2O masers trace the 
blue and red- shifte d lobes of the OF1 outflow (see Fig. [9] and 
ILeurini et aI1l2009l) . If we take into account (i) the EHV OF1 
red-lobe in CO (2 - 1), (ii) the SO red- and blue-shifted emis- 
sion, and (iii) the water maser spots, then the contribution of OF1 
to the kinematics of the CH3CN (I24 - 11 4) transition appears 
reasonably evident. Furthermore, Fig. |9]also shows red-shifted 
emission elongated towards the R2 red lobe of OF2-OF3 and 
blue-shifted emission close to the B 1 peak of OF2-OF3 as ob- 
served in CO (see Fig. [TJ. In Fig. [9] we also show the velocity 
field of the CH3CN (I23 - II3) line which has been excluded 
from our analysis because of weak blending with other features 
towards the centre position. However, being lower in energy than 
the (I24-I I4) transition, the (I23-I I3) line has a larger emitting 
size and can help the interpretation of the kinematics of the other 
CH3CN transitions. In this case, the N-W velocity gradient seen 
in the high J CH3CN transitions and to some extent also in the 
(I24 - 1 14) line is less pronounced, while red-shifted emission is 
clearly detected towards the Rl position in OF1 and towards the 
red-shifted lobe of the OF2-OF3 outflows. 

As an additional check, we derived position-velocity plots 
along several slits placed to sample all possible directions: the 
OF1 axis (P.A. -50°), the NE-SW direction (P.A. 30°, as rep- 
resentative of the three possible directions of the OF2-OF3 out- 
flows), two slits with position angles 60° and -20°, the N-S and 
E-W directions. All six slits are shown in Fig. [9] In Fig. QT| 
we show the PV diagram along the OF1 axis which shows that 
CH3CN (I24 - 1 14) emission is spread over a broad range of ve- 
locities (~ 15 km s~'). We then extracted spectra of the CH3CN 
(I24 - II4) transition for each slit, and fitted the line with a 
Gaussian profile. For each slit, we used a linear fit to repro- 



duce the peak velocity as function of the offset position along 
the slit. Results are shown in Fig. [12] The largest velocity gradi- 
ents are detected along the E-W direction (~ 114 km s _1 pc" 1 ) 
and the OF1 jet (~ 99 km s _1 pc~'). The gradients along 
OF2-OF3, and the slits at 60° and -20° have similar gradients 
(32 - 55 km s pc~'). Only the velocity gradient along the N-S 
direction is negligible, and no fit could be performed for this slit. 

Clearly, the velocity gradient along OF1 plays an important 
role in the kinematics of the inner gas around the protostars . This 
strongly supports our interpretation that the velocity field of the 
CH 3 CN (12 4 -1 1 4 ) line is affected by the OF1 outflow. In Sect.|6] 
we will discuss different scenarios for the velocity gradient de- 
tected along the E-W axis. 

The velocity field of the OCS (19-18) transition (Fig.Qj) is 
consistent with those of CH3OH and of the CH3CN (12 4 - 1 1 4 ) 
transition and shows red-shifted emission towards the red peak 
Rl of OF1 and R2 of OF2-OF3 and blue-shifted emission close 
to one of the blue peaks of OF2-OF3, Bl. O b CS (not shown 
here) is consistent with the picture depicted by high J CH3CN 
transitions. 

Finally, the HNCO emission shows a velocity gradient in the 
NE-SW direction along the axis of OF2-OF3 (see Fig. fT4l>. The 
first moment distribution of the C2H5CN (not shown here) is 
consistent with those of the high J CH3CN lines. 

6. Discussion 



As seen in Sect. 15. 21 different transitions trace different kinemat- 
ics. On the base of their first moment maps we can distinguish 
three cases: 

1. transitions with a velocity gradient along the NE-SW direc- 
tion (OF2-OF3), and additionally red-shifted emission to- 
wards NW (OF1). This is seen in CH 3 OH, CH3CN (12 4 - 
11 4 ), and OCS; 

2. transitions with a velocity gradient along the NE-SW axis 
(OF2-OF3) (HNCO); 

3. transitions with a velocity gradient along the E-W axis (high 
excitation CH3CN lines). 

As seen before, the first case can be explained by a combi- 
nation of the blue- and red-shifted lobes of the OF1, OF2-OF3 
outflows, indicating forCH 3 OH, CH 3 CN(12 4 - 1 1 4 ), and OCS at 
least a non negligible contamination from the multiple outflows 
in the IRAS 17233-3606 cluster. 

While the association of CH3OH and OCS with molec- 
ular outflows is not surprising since they have been de- 
tected in several chemical active molecular outflows (e.g., 
iBachiller & Perez Gutierrez! 119971) . the detection of a simi- 
lar behaviour in the CH3CN (I24 - 11 4) transition is unex- 
pected. CH3CN is usually considered a pure hot core tracer 
and is often used to investigate velocity fields a round massive 
YSOs to search for evidence of r otatio n (see ICesaroni et al.l 
1 19991: iBeltran et all l2004al I2005L 1201 ll b ut also discussion 
in iGibb et all |2004| and | Arava et al.l 120081). Howeve r, in the 



last years lArce et al.l (120081) and ICodella et al. (2009) detected 
CH3CN towards the molecular outflow driven by the low-mass 
YSO Ll 157-mm, hence showing that also this species can be al- 
tered by the passage of bow shocks, although with abundances 
much lower than in hot cores and hot corinos around high- and 
low-mass YSOs (e.g^ , lOlmi et all 1 19931; IBeltran et all 120051; 
iBottinelli etail 120071) . Additionally. IZapata et all (l2010ah ~re- 
cently proposed that CH3CN emission in Orion-KL is also ex- 
cited by shocks. In the present case, the velocity gradient de- 
tected in CH3CN could be due to outflow-envelope interactions 
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a [2000] 

Fig. 9. Summary of the high angular resolution observations available for the IRAS 17233 cluster (see previous figures for details 
and references). Left: Maps of the H2 emission (grey and solid contour levels). The blue and red contours represent the SO(5e - 4s) 
blue- and red-shifted emission; the dashed red contour is the CO(2 - 1) red-shifted emission at extremely high velocities. Overlaid 
on the map are also the positions of the blue and red OH (squares) and H2O (triangles) masers. Middle: Overlay of the CH3CN(124- 
1 14) first moment map (colour scale) with the H2 emission map (black contours). The black lines represent the directions used in 
the position-velocity plots of Fig.[T2]and are: the OF1 axis; the OF2-OF3 axis (NE-SW line); the N-S and E-W directions; the cuts 
with position angles of -20 and 60°. Right: Overlay of the CH3CN(123-1 13) first moment map (colour scale) with the H2 emission 
map (black contours). The black lines represent the axes of the OF1 and OF2-OF3 outflows. The Rl, R2 and Bl positions in OF1 
and OF2-OF3 respectively are marked by black squares. 



as already seen in low-mass star class sourc es dBeltran et all 
2004b; lArce & Sargeiitl 12001 iLee et all 12001 in other molecu- 
lar species. 

The second kinematical pattern (HNCO) can be easily as- 
sociated with the motions of the OF2-OF3 outfl ows. HNCO is 
abundant in densest regions o f molecular clouds (I Jackson et all 
19X4: IWvrowskietaTII1999h as well as in regions of low ve- 
locity shocks dlvlartm et all 120081) . However, very recently also 
HNCO h as been detected towards the che mical active LI 157 
flow by [Rodriguez-Fernand ez et al.l (1201 Oh . showing that the 
HNCO line profiles exhibit the same characteristics as other 
well-known tracers of shocked gas and supporting the idea that 
HNCO is a good tracer of interstellar shocks. Our data support 
this scenario, reporting the first detection of HNCO emission in 
outflows driven by massive YSOs. 

The third case (velocity gradient almost perpendicular to 
OF2-OF3) is the most intriguing since CH3CN emission has 
been extensively used in the last years to trace rotating toroids 
supposed to host accretion disks around massive YSOs. This has 
been used to support a mechanism for the formation of massive 
stars similar to that of low-mass stars. In the case of IRAS 17233, 
at least three possible interpretations can be given for the E-W 
velocity gradient detected in the CH3CN lines: rotation, contam- 
ination by outflows, velocity differences between cores. 

In the first case, the E-W velocity gradient could trace a ro- 
tating structure perpendicular to OF2-OF3 , the dominant out- 
flow o f this region. In this case, following [Beltran et al. (2004a, 
120051) we computed the dynamical mass Md yn assuming equilib- 
rium between centrifugal and gravitational forces. We took as 
rotation velocity V ro , = 3 km s , value corresponding to half 
the velocity range measured in the CH3CN (H4 - 1 14) line and 
as radius half the deconvolved FWHM size of the continuum 



emission (S - 5'.'3, see Table |3). This corresponds to a dy- 
namical mass of 27 M for an edge-on toroid. This is likely a 
lower limit to the real dynamical mass since the inclination an- 
gle is unknown. Such value is consistent with those estimated by 
iBeltran et all d2004all201 ll) for other sources. 

However, contamination by outflows questions this interpre- 
tation, since non negligible velocity gradients are detected along 
the OF2-OF3 axis and along other directions (see Sect. I5.2l and 
Fig. [T2l . In this context, the velocity gradient along the E-W 
direction could be due to a composition of motions along OF2- 
OF3 and OF1. Finally, since four cm compact VLA-sources are 
in the beam of our SMA data, and they are roughly located 
along the E-W axis, one could consider that the kinematical pat- 
terns here reported are simply due to a velocity difference of the 
sources. 

Clearly, higher angular observations are needed to discrim- 
inate between these interpretations. Given its low declination, 
ALMA is the ideal instrument to investigate the source at higher 
resolution. 



7. Implications for studies of high-mass YSOs 

Although not conclusive for the search of tracers of rotation 
around massive YSOs, we believe that our study of IRAS 17233 
has important implications for searches of disks in high-mass 
star forming regions. IRAS 17233 is an exceptionally close and 
active site of massive star formation. Therefore, even observa- 
tions with a limited angular resolution allow us to resolve small 
linear scales. Since the distance of the source is D < 1 kpc, our 
beam of ~ 5.4" x 2" corresponds to a maximum linear scale 
of 5400 x 2000 AU, thus comparable with state-of-the-art in- 
terferometric observations of more distant massive YSOs. For 
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example, the recent study of Beltran et al] (1201 ll) was performed 
with an linear resolution between 2400 and 50 000 AU. 

To verify whether the detection of OF1 in the CH3CN 
(I24 - 1 14) transition would have been possible with lower lin- 
ear resolution (i.e., a more distant object), we smoothed our data 
cube of a factor of two in spatial resolution. This corresponds to 
observations of a massive YSO at a distance twice as large as 
that of IRAS 17233, but with the same angular resolution as our 
observations of IRAS 17233. We then computed moment maps 
for the CH3CN (I24 - II4) transition following the same pro- 
cedure described in Sect. [5] We also smoothed the CO (2-1) 
data cube to the same angular resolution, and derived integrated 
intensity maps for the blue- and red-shifted emission in the same 
velocity range as in Fig. Q] The resulting maps are presented in 

Rg. [B] 

The data at lower angular resolution show a compact re- 
gion detected in CH3CN with a velocity gradient perpendicu- 
lar to the dominant outflow of the region, OF2-OF3. Clearly, 
the red-shifted emission along the OF1 axis is not longer de- 
tected, and the derived velocity field is very similar to that de- 
tected for the CH3CN {K = 7, 8) lines. In other words the picture 
of Fig. [l5]could easily be explained with a pure velocity gradi- 
ent perpendicular to the N-S outflow as seen in CO. We note 
that for sources at distances > 10 kpc also the detection of H2 
is challenging. Indeed, in the case of IRAS 17233, the H2 data 
were fundamental to reveal OF1, since H2 gives a much better 
picture of the outflow multiplicity than CO because it traces hot 
jets and is not contaminated by infalling envelopes and/or swept 
up cavities. 

8. Conclusions 

In this study, we analysed interferometric line observations of 
the massive star forming region IRAS 17233 at 230 GHz with an 
angular resolution of 574 x 179. We studied the spatial distribu- 
tion and velocity field of different lines and compared them to the 
distribution of outflows in the region to investigate the nature of 
the emission. We report the first detection of HNCO in molecu- 
lar outflows from massive YSOs. Although other interpretations 
cannot be ruled out, we find that the velocity gradient observed 
in CH3CN could be significantly affected by the OF1 outflow. 
This follows the recent detection of both molecules in the chem- 
ically active outflow Ll 157, powered by a low-mass YSO. The 
possible detection of CH3CN in molecular outflows questions 
the choice of the tracer to probe rotating structures around YSOs: 
one has to keep in mind that the same molecule could be trac- 
ing outflows in some objects and rotating disks/toroids in others. 
Further studies of high-mass star forming regions are needed to 
verify whether IRAS 17233-3606 represents a standard labora- 
tory where to investigate star formation or whether it is a spe- 
cial case with a too complex scenario. The present observations, 
with the possible association of CH3CN with the OF1 outflow 
in a few thousands AU around the peak of the continuum, stress 
the importance of (i) very high-spatial resolution, less than 1000 
AU, and of (ii) H2 images for detailed studies of kinematics in 
regions as those around massive YSOs. Assuming typical dis- 
tances of massive star forming regions between 5 and 10 kpc, a 
resolution of 1000 AU corresponds to an angular size of 071- 
072. Even smaller angular resolutions are then needed to study 
in details the kinematics of such regions. Only studies of mas- 
sive YSOs with angular resolution < 071 (for distances > 5 kpc) 
will clarify whether or not IRAS 17233 represents a special case 
and will shed light on to the choice of the molecular tracers for 
rotating structures. 



As a final remark, given the chemical differences in hot 
cores, and contamination by outflow emission, potentially in all 
species, it well may be the case that no molecular species or 
transition will work for every source under every circumstance. 
However, we can try to formulate criteria for convincing detec- 
tion of a rotating disk/toroid in millimeter molecular emission: 
(i) such a detection should be based on the analysis of high ex- 
citation multiple lines and multiple species; (ii) the lines used in 
the analysis should be unblended; (iii) any possible contribution 
from outflows should be ruled out; (iv) all lines should be con- 
sistent with a single velocity field of a rotating body. As noted in 
the previous section, selective tracers of hot jets as H2, or SiO, 
are better suited to understand the multiplicity and morphology 
of outflows than CO and may therefore better help to understand 
the kinematics of the cores. 
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Fig. 4. Integrated emission of selected transitions of CH3OH in 
the region above the 5cr detection in the line data cubes (~ 0.5 
Jy beam -1 ). The black solid contours show the 1.3 mm contin- 
uum emission as in Fig. [2] The solid line shows the axis of OF1, 
while the dashed line has been chosen, for sake of clarity, to rep- 
resent the three possible directions of the N-S outflows drawn 
in Fig.Q] The red squares mark the Rl, R2 and Bl peaks of the 
CO (2 - 1) emission detected towards OF1 and OF2-OF3 (see 
Fig. |l(a)| >. The white ellipse represents the FWHM 2-D Gaussian 
fit to the data. 
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Fig. 5. Integrated emission of the OCS (19-18) line in the region 
above the 5<x detection in the line data cubes (~ 0.5 Jy beam" 1 ). 
The black solid contours show the 1.3 mm continuum emission 
as in Fig. |2 Symbols are as in Fig. [4] 
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Fig. 6. Integrated emission of selected transitions of CH3CN in 
the region above the 5<x detection in the line data cubes (~ 0.5 
Jy beam" 1 ). The black solid contours show the 1.3 mm contin- 
uum emission. Symbols are as in Fig. |4] 
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Fig. 7. Integrated emission of selected transitions of HNCO in 
the region above the 5<x detection in the line data cubes (~ 0.5 
Jy beam 1 ). The black solid contours show the 1.3 mm contin- 
uum emission. Symbols are as in Fig. [4] 
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Fig. 8. First moment maps of selected transitions of CH3OH. 
The black solid contours show the 1.3 mm continuum emission. 
Black squares mark the Bl, Rl and R2 positions of Fig. |l(a)| 
The other symbols are as in Fig. |4] The systemic velocity of the 
source lies between -3.4 and -2.5 km s _1 (see Sect. 15.2b . 
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Fig. 10. First moment map of the CH 3 CN (12 7 - 1 1 7 ) line. The 
black solid contours show the 1.3 mm continuum emission. 
Black squares mark the Bl, Rl and R2 positions of Fig. |l(a)| 
The other symbols are as in Fig. |4] The systemic velocity of the 
source lies between -3.4 and -2.5 km s _1 (see Sect.P 
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Fig. 12. Peak velocity of the of CH 3 CN(12 4 -1 1 4 ) line along six 
different directions outlined in Fig. [9] as function of the offset 
position along the slit. Different colours denote different slits. 
The solid lines represent the fit to the data. 
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Fig. 11. Position-velocity plot of the CH 3 CN(12 4 - 11 4 ) 
line along the OFl axis. Contour levels range from 5cr 
(0.5 Jy beam -1 ) in step of 20<x. The horizontal line marks the 
systemic velocity of -3.4 km s , the vertical line the crossing 
point of all slits. 



Fig. 13. First moment map of the OCS (19 - 18) transition. 
The black solid contours show the 1.3 mm continuum emission. 
Black squares mark the Bl, Rl and R2 positions of Fig. |l(a)| 
The other symbols are as in Fig. |4] The systemic velocity of the 
source lies between -3.4 and -2.5 km s~ ] (see Sect. 15. 2b . 



13 



S. Leurini et al.: On the kinematics of massive star forming regions: the case of IRAS 17233-3606 




[2000] 



Fig. 14. First moment maps of the HNCO (10 - 9 ) line. The 
black solid contours show the 1.3 mm continuum emission. 
Black squares mark the Bl, Rl and R2 positions of Fig. |l(a)| 
The other symbols are as in Fig. |4] The systemic velocity of the 
source lies between -3.4 and -2.5 km s" 1 (see Sect. 15.21 ). 
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Fig. 15. Velocity gradient of the CH 3 CN (12 4 - 11 4 ) transition 
smoothed to an angular resolution of 10786 x 3778. The red 
and blue solid contours show the red- and blue-shifted CO(2- 
1) emission smoothed to the same angular resolution. 
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